Dynamical Correlations of the Kagome 5 = 1/2 Heisenberg Quantum Antiferromagnet 



Andreas M. Lauchli 1 and Claire Lhuillier 2 

' Max Planck Institut fiir Physik komplexer Systeme, Nothnitzerstr. 38, D-01187 Dresden, Germany. 
2 Laboratoire de Physique The'orique de la Matiere Condensee, UMR 7600 CNRS, 
Universite Pierre-et-Marie-Curie, Paris 6, 75252 Paris cedex 05, France. 
(Dated: January 8, 2009) 

We determine dynamical response functions of the S = 1/2 Heisenberg quantum antiferromagnet on the 
kagome lattice based on large-scale exact diagonalizations combined with a continued fraction technique. The 
dynamical spin structure factor has important spectral weight predominantly along the boundary of the extended 
Brillouin zone and energy scans reveal broad response extending over a range of 2 ~ 3J concomitant with 
pronounced intensity at lowest available energies. Dispersive features are largely absent. Dynamical singlet 
correlations - which are relevant for inelastic light probes - reveal a similar broad response, with a high intensity 
at low frequencies u/J < 0.2 J. These low energy singlet excitations do however not seem to favor a specific 
valence bond crystal, but instead spread over many symmetry allowed eigenstates. 

PACS numbers: 75.10Jm, 75.40.Mg, 75.40.Gb 



Introduction The 5=1/2 Heisenberg quantum antifer- 
romagnet (AFM) on the kagome lattice is a key model for our 
theoretical understanding of highly frustrated quantum mag- 
nets in two spatial dimensions. In contrast to other model sys- 
tems, such as the checkerboard magnet, it continues to hide 
the true nature of its ground state despite a longstanding ef- 
fort. On the experimental side many kagome like materials 
were discovered and characterized over the years, and some 
of them seem to get close to the goal of an experimental real- 
ization of a perfect Heisenberg 5=1/2 kagome system. Re- 
cently the synthesis of powder samples of the Herbertsmithite 
ZnCu3(OH)eCl2 [1] and the subsequent experimental investi- 
gations sparked a new wave of theoretical interest in this long 
standing problem [2]. 

So far most of the theoretical studies focused on ground 
state properties, discussing various possible phases, such as 
valence bond solids [3], gapped spin liquids of different 
kinds [4, 5, 6], and also stable critical phases [7, 8]. Much 
less attention however has been paid to the precise nature 
and form of the low-energy excitations visible in frequency 
resolved probes such as inelastic neutron scattering or light 
scattering techniques. In the present Letter we fill this void 
and present a detailed numerical study of the dynamical re- 
sponse of 5 = 1/2 kagome systems in the spin triplet and 
singlet channels and formulate predictions to be tested in in- 
elastic scattering experiments. 

We study the 5=1/2 Heisenberg quantum antiferromag- 
net on the kagome lattice, governed by the Hamiltonian: 

H = J^S i -S j (1) 

where J > is the antiferromagnetic nearest neighbor ex- 
change coupling. Our results are based on large-scale exact 
diagonalizations of up to N = 36 spins, supplemented by 
the continued fraction method [9] for dynamical correlations 
functions. For the dynamical quantities we performed 500 
(3000) iterations in Hilbert spaces of dimensions 3.8 x 10 s 
(4.5 x 10 9 ) for the spin and singlet dynamics respectively. 



A lot of our present understanding of 5 = 1/2 systems is 
based on a series of exact diagonalization studies [10], which 
convincingly showed the absence of magnetic order, and re- 
vealed a puzzlingly high density of low-energy singlet and 
triplet excitations, a so far unique phenomenon. In the follow- 
ing we explore how the huge density of low-lying excitations 
affects experimentally relevant response functions. 

Static spin response The static spin structure factor is 
given by the following expression: 

5*(Q) = i£e-^SJ, 

j 

5(Q) = (5 2 (-Q)5 z (Q)) , (2) 

where the wave vector Q is not restricted to lie in the first 
Brillouin zone (BZ). 

The numerically determined static spin structure factor for 
the 36 sites sample is shown in Fig. 1, as an intensity plot 
covering multiple BZs [right hand side, plot (1)], and along 
a path in the extended BZ [plot (2)]. The response is strong 
and broad along the zone boundary of the extended BZ, re- 
vealing the short ranged nature of the antiferromagnetic cor- 
relations [10]. In the N = 36 groundstate there are additional 
small peaks at the point (g), a feature which is also reported in 
a recent DMRG study [6]. Comparison with a N = 24 sample 
(not shown) confirms that the broad response along the bound- 
ary of the exended BZ is a generic, size-independent feature, 
characterizing a state with spin correlations which are decay- 
ing rapidly beyond the nearest neighbor sites. To contrast this 
result with a magnetically ordered state we show the static re- 
sponse in the q = [Fig. 1(3)] and ^3 x ^ [Fig. 1(4)] states, 
which have been obtained using a strong J2 coupling with the 
appropriate sign. 

Dynamical spin structure factor The energy and momen- 
tum dependence of the dynamical structure factor: 

5(Q,w) = --Im(5*(-Q) — 1 5 Z (Q)) , 

7T UJ - (H - E GS ) + IJ] 

(3) 
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FIG. 1: (Color) Dynamical spin structure factor of the N=36 sample. The eight panels display frequency scans 5(Q, u>) (g = 0.02 J) at 
labeled wavevectors Q in the extended Brillouin zone shown in the lower right center. Note that the intensity scales differ among the different 
panels. The T point has no weight and is not shown. The blue vertical lines show the pole location and intensity of the continued fraction. The 
vertical dotted magenta line denotes the finite size spin gap in the corresponding momentum sector. The dashed red line marks the position 
of the first frequency moment uj — J du lo 5(Q, lo) /<S(Q). In the rightmost column the static spin structure factor of the pure Heisenberg 
model on the kagome lattice is shown, as an intensity plot (1) and along the path V — (e) — (g) — T (2). The static structure factor for the 
q = (3) and \/3 x \/3 (4) Neel order states induced by appropriate second neighbor couplings are also displayed. 



is directly relevant for inelastic neutron scattering (INS) ex- 
periments and therefore a quantity of central interest. In mag- 
netically ordered systems we expect to see dispersive, long- 
lived spin waves [11], while one-dimensional systems in ap- 
propriate regimes can display spinon continua with a rich 
structure [12]. 

Our numerical results for the N — 36 kagome lattice are 
presented in the left part of Fig. 1 . The shaded panels dis- 
play an energy scan at the wave vector indicated by the panel 
position and its label, referring to specific points in the ex- 
tended BZ. Each panel displays the broadened (77=0. 02J) 
spectral function (black line), the locations and weights of 
the poles of the continued fraction expansion (blue vertical 
lines), the finite size spin gap in the corresponding momen- 
tum sector (dotted vertical line), and the first frequency mo- 
ment ui(Q) — J dwwiS(Q,w)/5(Q) (dashed vertical line). 

Consistent with the static structure factor [by virtue of the 
sum rule S(Q) = J ' dui S(Q,ui)], the dynamical spin re- 
sponse function concentrates essentially along the boundary 
of the extended BZ. The main feature of this system is the 
stretching of the magnetic response in each Q-sector on a very 
large number of excited states spanning a large bandwidth 



of 2 ~ 3 J, starting immediately above the (finite-size) gap. 
Furthermore there seems to be a pronounced enhancement 
of the intensity at small u. The different spectral functions 
look rather similar, suggesting an approximate factorization 
S(Q, lo) ~ S(Q) x f(uj), at least at intermediate and high lo. 
The overall picture is definitely quite different from the spec- 
trum of a Neel ordered system on the same system size, where 
an overwhelming part of the spectral weight is carried by very 
few poles in each Q-sector associated to the Bragg peak and 
the one-magnon modes respectively [13]. Still at some wave 
vectors the lowest pole carries significant weight [especially 
at (g)]. The origin of this feature remains to be elucidated but 
could potentially come from an algebraic divergence in one 
scenario [8] or from triplon excitations on top of (remnants 
of) a valence bond crystal in a different scenario [14]. 

In order to address finite-size effects we present two spec- 
tral functions at the wave vectors (g) and (i) for N = 24 
and 36 spins in Fig. 2(a) and (b). The characteristic width in 
energy as well as the prominent response at low lo for wave 
vector (g) are clearly stable with respect to finite size ef- 
fects. In panels (c) the local dynamical spin correlation func- 
tion S\ oc (lj) oc J dQ S(Q,lo) is shown for different system 
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FIG. 2: (Color online) Upper panel: Finite size behavior of the 
spectral functions at two different points in the Brillouin zone: g (a) 
and i (b) in the convention of Fig. 1. Lower panel: (c) Local spin 
autocorrelation function Si oc (uj) oc J dQ 5(Q, uj) for N = 24, 30 
and 36 sites, (d) The cumulative spectral weight as a function of 
Lu/ J. All spectral functions have been broadened using r\ — 0.05 J. 



sizes, again highlighting the stability of the overall shape of 
the spectral function. Finally panel (d) presents the cumula- 
tive spectral weight as a function of uj/ J, revealing e.g. that 
uj ~ 0.95 J is the median frequency for the total spin response. 

Effect of impurities We have studied the influence of a low 
concentration of nonmagnetic impurities on the spin dynam- 
ics by depleting a N = 27 sample by one site. The averaged 
dynamical spin response closely resembles Fig 2(c) with an 
additional resonance-like feature at uj ~ J due to the strong 
singlet forming on the two bonds next to the vacant site [15]. 
Due to its local nature, this feature is expected to be geneti- 
cally present in 6>(Q, uj) as well. 

Singlet fluctuations In order to assess the importance of 
the abundant number of low energy singlet excitations for op- 
tical probes and to investigate the tendency towards valence 
bond crystal ordering, we study the local dynamical fluctua- 
tions of a nearest neighbor dimer operator: 
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The interest in this quantity is twofold. First Eq. (4) is closely 
related to the Raman or RJXS reponse of a spin system and 



FIG. 3: (Color) Dynamical singlet fluctuations [Eq. (4)] for three 
different systems. Main plot: kagome lattice. Upper inset: Checker- 
board lattice with a plaquette-like valence bond crystal ground state. 
Lower inset: Unfrustrated square lattice exhibiting Neel order. The 
plotted quantity represents well the qualitative features of the Raman 
response of the three systems. 



thus reveals the qualitative features of the inelastic light scat- 
tering response. And second, we expect a spontaneous trans- 
lational symmetry breaking due to dimerization to show up as 
an important ui — > contribution. 

The fluctuation spectrum for the N = 36 kagome system 
is shown in Fig. 3, where a broad response from the lowest 
singlet up to energies ~ 4J is seen together with a strong 
increase of the response towards the lowest energies. The 
kagome result can be compared to the response of the Heisen- 
berg model on a checkerboard lattice (upper inset in Fig. 3) 
and the square lattice (lower inset in Fig. 3), where in both 
cases the physical origin of the response is essentially under- 
stood. On the Neel ordered square lattice the dimer-dimer re- 
sponse is assigned to the two-magnon continuum with a maxi- 
mum strength around 3 J. This is obviously very similar to the 
Raman response of the square lattice Heisenberg AFM [16]. 
On the checkerboard lattice with its plaquette valence bond 
crystal ground state [17] the dimer-dimer response function 
shows different frequency domains: a single low-lying peak 
(shaded in red) originating from the valence-bond symmetry 
breaking partner of the ground state, followed by two domains 
at non zero-frequency, first a range of singlet excited levels 
(shaded orange) which have been convincingly explained as 
valence bond crystal domain-wall excitations [17, 18] and a 
second range uj > J corresponding to a multi-triplon con- 
tinuum [17] (shaded in brown). It is clear from this com- 
parison that the kagome lattice does not show typical valence 
bond crystal characteristics exemplified by the checkerboard 
magnet. Still the response on several levels up to uj < 0.2 J 
seems to be particular strong. However this response spreads 
on many low lying levels in any symmetry sector which can 
be excited by the dimer-dimer operator. We do not see any 
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clear precursor of a specific spatial symmetry breaking pattern 
when comparing the excited levels to different valence bond 
crystal symmetry predictions summarized in Ref. [19]. Possi- 
ble reasons are: i) the absence of VBC order (see Ref. [6] for 
a similar conclusion) or ii) a very weak ordering with a very 
large unit cell which has not still emerged from competing 
orders. 

Comparison to theoretical proposals Both the spin and 
the dimer dynamical fluctuations of a S = 1/2 kagome sys- 
tem are intrinsically broad, and are not easy to reconcile with 
the excitation spectrum scenarios for the various proposed 
ground states. For example the dynamical spin correlation 
functions lack the characteristic coherent triplon excitations, 
which are expected on top of a valence bond crystal [14]. 
On the other hand the spectrum is much too dense to be ex- 
plained by a spinon continuum with a energy scale of order 
J [4]. A critical spin liquid would possibly have similar low- 
oj enhanced response [8] as we report, but the current theoret- 
ical predictions would still predict dispersive structures with a 
bandwidth significantly larger than our numerical bandwidth. 
Interestingly our fully quantum mechanical results show some 
similarities with fluctuation spectra of classical highly frus- 
trated systems [20], where the macroscopic ground state de- 
generacy plays an important role. It is thus possible that the 
intermediate and high energy response of the T — kagome 
S = 1/2 quantum magnet resembles that of a classical co- 
operative paramagnet at finite temperature, while the precise 
low-energy response will ultimately be dictated by the yet un- 
known true nature of the ground state. 

Comparison to experiments Our results provide definite 
predictions for the spectral response of a perfect S = 1/2 
kagome system that can be checked in experiments. First of 
all a pure Heisenberg model on the kagome lattice has a static 
spin structure factor with a negligible weight inside the first 
BZ and most of the weight is concentrated on the boundary 
of the extended BZ. While neutron diffraction data on many 
S > 1/2 systems is in good agreement with this expectation, 
the same data for the 5=1/2 Herbertsmithite [21, 22] is diffi- 
cult to reconcile with a simple Heisenberg model and possibly 
points to some remnant magnetism that may be due to impu- 
rities or Dyaloshinsky-Moryia interactions. Similarly for the 
inelastic response: while S = 5/2 Deuterium-Jarosites are in 
nice qualitative agreement [23] with our results, the reported 
inelastic neutron signals for the Herbertsmithites [21, 22] is 
concentrated at energies lower than what we would expect for 
a pure Heisenberg model. Therefore more experimental and 
theoretical work is needed in order to understand the Herbert- 
smithite response. Finally our results might also provide a 
first step towards an understanding of the anomalous dynam- 
ical spin fluctuations revealed by /isr on SrCr 8 Ga 4 Oi9 [24] 
and related compounds. 

Perspectives Based on dynamical correlations functions 
in the tripled and singlet channels we have shown that the 
kagome system is a highly fluctuating magnet with broad re- 
sponse peaking at low energies. This fluctuating background 
provides also a natural explanation for the reported absence 



of quasiparticles in a doped kagome system [25]. Building 
on these results it will be interesting to understand the evo- 
lution of the dynamical response at finite temperature as well 
as the effect of a magnetic field. Given that many experimen- 
tal kagome systems have some (small) Dzyaloshinsky-Moriya 
interactions their influence on the dynamical response is also 
worth studying. 
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